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Abstract We propose a method for deterministic im- 
plantation of single atoms into solids which relies on a 
linear ion trap as an ion source. Our approach allows a 
deterministic control of the number of implanted atoms 
and a spatial resolution of less than 1 nm. Furthermore, 
the method is expected to work for almost all chemical 
elements. The deterministic implantation of single phos- 
phor or nitrogen atoms is interesting for the fabrication 
of scalable solid state quantum computers, in particular 
for silicon and diamond based schemes. A wide range 
of further applications is expected for the fabrication of 
nano and sub-nano electronic devices. 

PACS 03.67.-a; 29.25.Ni; 61.72.Ji; 81.16.Rf; 85.40.Ry 
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I> 1 Introduction 



A future solid state quantum computer (QC) may be 
based on implanted single ions, each of them carrying a 
single unity of quantum information (qubit). Logic gate 
operations rely on the well controlled coherent interac- 
tion between individual qubits. Any coupling to the solid 
state bulk material would destroy the -in general- entan- 
gled global quantum state of the device. This reasoning 
has led to the concept of qubits logic states which are 
encoded into hyperfine states of single atomic phospho- 
rous ions '^^P, embedded in a pure Si mono crystal 
More recently, experiments with single nitrogen- vacancy 
(NV) centers in diamond crystals have demonstrated sin- 
gle qubit operations T and two-qubit quantum logic 
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gate operations mediated by hyperfine coupling [3], and 
a long coherence time U as necessary for a QC according 
to the DiVincenco criteria [S]- 

However, the positioning accuracy of NV-defects was 
so far limited to 50 nm using a damage-implantation 
technique ■ The future quest is to fabricate an array 
of NV-centers with nm-precision for a scalable solid state 
QC. So far, no scheme for deterministic singly-charged 
single ion implantation has been yet reported for this 
spatial resolution. 

The paper is organized as follows: After a short dis- 
cussion of current approaches of ion implantation, we de- 
scribe the novel method for single ion implantation with 
sub-nm spatial resolution. This includes a Paul trap as 
an ion source and a micro Einzel lens as an ion opti- 
cal element for focusing through a pierced atomic force 
microscope (AFM) tip. Finally, we sketch the basic con- 
cepts of a solid state QC based upon NV color centers. 
Though this implantation method may be applied for a 
large variety of atomic species we will focus here on the 
specific case of N"*" ions being implanted into diamond 
for the deterministic generation of NV color centers. 



2 Ion implantation 

The nitrogen-vacancy center in diamond is an atomic 
point-defect with a structure that consists of a nitrogen 
atom and a vacancy in the next lattice position (Fig.^). 
It was demonstrated that NV defects can be created in 
nitrogen-rich (type lb) diamond by the creation of va- 
cancies and a subsequent annealing step 6 . During this 
annealing process, vacancies migrate towards nitrogen 
atoms forming NV defects. The major drawback of this 
approach is its poor positioning accuracy of defects due 
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Fig. 1 a) Atomic structure of the nitrogen-vacancy center 
in diamond consisting of a nitrogen atom (N) and a vacancy 
(V) in the next lattice position, b) Scheme of the electronic 
and spin energy levels. Excitation and fluorescence is in the 
visible at 637 nm. Ground state levels may serve to store 
quantum information. 




to the diffusion of vacancies (50 nm). Alternatively, NV 
defects can be created using nitrogen implantation into 
nitrogen-free diamond (type Ila) jHHH] ■ Ion beam implan- 
tation has the advantage that the preselected position of 
the N-atom now defines the place of the NV center. Only 
the vacancies migrate during annealing (T^^ 850°C) as 
they become mobile above 650° C. The diffusion of sub- 
stitutional nitrogen starts only above 1500°C. The cre- 
ation of a single NV center relies on the migration of a va- 
cancy to the nitrogen atom during the annealing process. 
The efficiency depends on the energy of the implanted 
ion and increases with higher energy as an increasing 
amount of vacancies is generated by the implantation 
process itself. 

As an illustration of the latter method we show an 
array of NV centers generated using focused N+ ions 
(2 MeV) implanted into type Ila diamond (see Fig. [^J. 
The efficiency of NV center creation at this energy is 
about 50% |H1 ■ The spot size is about 300 nm on the tar- 
get surface. Due to the stochastic nature of the ion beam 
some spots contain more than one NV center, while oth- 
ers contain none. Additionally, because of the high ki- 
netic energies of the ions the observed spot size is mostly 
limited by ion straggling in the diamond matrix. 

However, for single ion implantation with an aspired 
resolution of less than 2 nm it is necessary to reduce the 
ion energy down to 1 keV in order to avoid lateral dis- 
placement caused by straggling^. Recent results indicate 
that at these low implantation energies the efficiency of 
the creation of a NV defect during the annealing process 
is decreased to less than 10%^. Subsequent implanta- 
tion of carbon ions could be used to increase the local 
density of vacancies near the nitrogen atom. Also other 
elements could be used which immediately generate color 
centers without the need of an additional vacancy like 
boron in diamond or chromium in sapphire. Regardless 
of the used ion species, the alignment of the ion beam is 



Fig. 2 Confocal fluorescence image of implanted NV centers 
generated from ions implanted into type Ila diamond. 
Color centers are excited with light near 637 nm and the 
emitted fluorescence is detected by using a confocal micro- 
scope (intensity encoded in color scale). 
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Fig. 3 Schematic of the setup used for deterministic implan- 
tation of single ions. The ions are trapped and subsequently 
laser cooled inside a Paul trap. A voltage change at the ex- 
traction electrode accelerates the ion towards the substrate. 
The ions are focused onto the substrate with an Einzel-lens 
mounted together with an AFM-tip used for alignment. 



realized by a pierced AFM-tip acting as a stencil nano 
mask. The use of a pierced AFM-tip as a collimator is 
already successfully established [111112) . However, this 
technique has two disadvantages: First, a pinhole size 
of 5 nm and below will produce proximity effects, affect- 
ing the spatial resolution. Secondly, this technique needs 
highly charged ions for the detection of the impact of 
one single ion. But highly charged ions may cause severe 
damage to the surface |13j , possibly destroying the solid 
state QC coherence. 

A more suitable method presented in this paper is 
the use of an ion trap as a single ion point source. This 
technique avoids any proximity effects and does not re- 
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quire any additional detection system. A pierced AFM 
tip acts as alignment tool for the ion lens with respect 
to the target. 

The basic ideas of the proposed system are shown in 
Figure 13 A linear Paul trap is used as a source of single 
ions with supreme performance from an ion-optical point 
of view. Single ions are trapped, cooled and extracted 
with a spatial- and momentum-uncertainty near to the 
fundamental limit solely given by quantum physics. Ad- 
ditionally, the resonance fluorescence of those ions can 
be routinely excited by laser radiation and observed with 
a camera system. Thus, the number of trapped ions is 
well known and a deterministic implantation of a given 
number of atoms is possible. The ions are extracted and 
focused down to the surface of the substrate with sub-nm 
resolution using a micro Einzel lens. Advantageously, it 
is combined with a pierced AFM tip for the exact place- 
ment and accurate positioning of the system. In order 
to realize this design the following questions have to be 
addressed: 

— Is it possible to extract one ion from a Paul trap 
without losing the position information? 

— What kind of ion optics is adequate to focus ions 
down to nm resolution? 

— How to manage the adjustment procedure for single 
ion implantation? 

— What is the expected final spatial resolution of the 
entire system? 



3 A linear Paul trap as a deterministic single 
ion source 

Paul traps are widely used for single ion experiments 
|14lll5lll6j . In the specific case of a linear Paul trap, the 
radio-frequency (RF) voltage is applied to a set of radial 
electrodes and generates a pseudo potential of several 
eV depth. The axial confinement is provided by addi- 
tional DC electrodes. Various forms and sizes of linear 
Paul traps are known: Typically, the RF-electrodes are 
in the shape of four rods [13 or four blades jTHj. The 
applied voltages result in a harmonic trapping poten- 
tial for all spatial directions. Trap sizes range from a 
few mm to a few 100 /im in case of small linear traps 
|19j . Correspondingly, the vibrational trap frequencies 
vary from several MHz down to 100 kHz. Due to laser 
Doppler cooling, the kinetic energy of the trapped ions 
can be reduced down to a few mK (corresponding to 
lO"'^ eV). The typical wave packet extension is 20 nm 
for a "'^'*N+ ion in the ground state at a trap frequency of 
Wtrap/27r = 2 MHz. Routinely, the ground state of vibra- 
tion is reached by more sophisticated laser cooling tech- 
niques such as resolved side band cooling on quadrupole 
transitions [20] or ground state cooling exploiting elec- 
tromagnetically induced transparency (FIT) [21]. As a 
result, the position- momentum uncertainty of the ion 
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Fig. 4 String of 10 *°Ca.+ ions plus 3 impurity ions excited 
and cooled by laser radiation near 397 nm and 866 nm. The 
ion crystal's fluorescence near 397 nm is imaged onto a CCD 
system. In the top image one impurity ion is in the center, the 
other two are at the far right. The average inter-ion spacing 
is 10 /im (axial oscillation frequency uJz/2tt— 108 kHz, radial 
frequency curad/^n = 1.4 MHz). Exposure time is 1 s. Between 
the upper, middle and lower image, the ion crystal was melted 
by tuning the laser radiation near 397 nm to the blue side of 
the atomic resonance. The positions of the impurity ions are 
changed after cooling the ion crystal again. Figure from |26| 
with kind permission of the author. 



wave function is limited only by the Heisenberg uncer- 
tainty principle. 

However, ion species which might be interesting for 
doping applications can usually not be laser cooled. This 
is due to the fact that the corresponding transition wave- 
lengths are extremely hard to achieve by laser sources or 
the level scheme does not feature an adequate two-level 
system suitable for an optical cycling transition. Unfor- 
tunately, this is the case for 3ip+ and for any isotope of 
nitrogen. This difficulty may be circumvented by sym- 
pathetic laser cooling [^28112411^ . We propose to sym- 
pathetically cool either phosphorous or nitrogen doping 
ions inside a string of ^''Ca"'" ions which are laser cooled. 
The linear crystal of ''^"Ca"'" ions will be observed with a 
CCD-camera whereas the doping ions show up as dark 
sites |26II23| (Fig. 21 . After Doppler cooling, a second 
cooling step via FIT will allow to reach almost the vi- 
brational ground state of all common modes of vibration 

E3 

Then, a single dark doping ion is separated together 
with a single '*'^Ca"^ ion from the entire crystal in order 
to determine the mass of the doping ion by measuring 
the vibrational frequency of the mixed crystal [221 • The 
observation of dark sites in the fluorescence, the ground 
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state cooling, and the deterministic separation of the 
doping ion are the key-elements of our method. 

Finally, the ion is extracted for implantation, being 
delivered with an ideal brightness at the physical limit. 

4 Point source ion optics 

The ion trap for single ions defines a completely new 
starting point for ion beam optics offering an ideal point 
source perfectly suitable for deterministic ion ray trac- 
ing. The theoretical accuracy limit of the point source is 
given by the uncertainty of the initial velocity distribu- 
tion after optical cooling. Beam broadening by stochastic 
space charge effects is completely avoided. 

For the numerical simulation of the ion's trajectory 
we first model the quasi-electrostatic potential of the 
Paul trap and then simulate the trajectories by ray trac- 
ing. In order to cross-check the results, two fully inde- 
pendent software packages have been developed. 

Method A relies on the calculation of the potential 
distribution of the trap and extractor lens using a finite- 
differential method (FDM) combined with a multi-grid 
technique to find starting values for a simultaneous over- 
relaxation (SOR) method |^. The calculation of the 
trap potential is performed on a 64x64x256 grid. The 
grid distance is chosen to be 100 fim. A ray tracing pro- 
cedure uses the resulting field distribution to calculate 
the electric force acting on the ion using a 5th order 
Runge-Kutta formalism with adaptive step size control 
PS) . A Monte-Carlo technique simulates the influence of 
starting parameters. The temperature is set as starting 
parameter for the random choice of the initial ion veloc- 
ity. Method B applies the commercial program Femlab^ 
for calculating the 3D trap potentials in an adaptive non- 
uniform grid. The trajectories are simulated in C-I-+ in 
a 8th order Runge-Kutta Prince-Dormand method^. 

In both simulations, the type and size of the linear 
quadruple ion trap is adapted from pMTTj . The trap 
consists of four parallel rods with 2 mm spacing and a 
length of 10 mm. Diagonally opposed rods are connected 
to a RF source with Lo^jl-K of 18 MHz and V^p = IkV, 
generating a quadrupole RF field. The axial confinement 
is provided by two disc shaped end caps at a distance of 
10 mm (outer diameter 6 mm) with a centered extraction 
hole of 1 mm diameter. The DC voltage at the end caps 
is set to 500 V. The parameters for the end caps are 
slightly modified to [T7] . 

To use the trap as an ion source all potentials are 
biased by -1-500 V and the temperature of the ions is as- 
sumed to be around 0.1 mK. The extraction of the ion is 
initiated by switching one of the end caps from 1000 V to 
V. Both simulation methods show that sub-nm resolu- 
tion (<1 nm) can be achieved with laser cooled ions. In 
the following we will focus only on the results of method 

^ Femlab; http;//www. comsol.com/ 

^ GSL: , http://www.gnu.org/software/gsT7] 



B. The rays of the ions inside the trap after extraction 
are shown in Figure[Sl One of the important questions for 
the optimization is the handling of the RF-ficld during 
the extraction procedure. Our first guess was to contin- 
uously apply the RF-field during extraction in order to 
maintain the radial confinement (Fig. [S] dashed). The 
lateral spread of the ion trajectories at a distance of 
10 mm from the trap center in the extraction direction is 
about ±50 nm. In order to focus the trajectories down 
to one spot at the target we simulate a small electro- 
static Einzel lens between trap and target, located at a 
distance of 90 mm from the trap center. The symmetric 
Einzel lens consists of a stack of three equidistant con- 
ducting layers separated by about 150 \xm. The middle 
layer is biased to -1-1000 V whereas the other layers are 
grounded. The resulting focal length for 500 eV ions is 
about 300 /im (Fig.O dashed). The spread of the impact 
positions of the ions at the target plane (Fig.jTla) is about 
20 nm which would be still one order of magnitude too 
high. If the position of the target plane is changed within 
1 /im the spot size hardly changes (Fig.|7^-c). This situ- 
ation cannot be improved by further lowering the initial 
temperature of the ions in the trap. The spread is caused 
by axial velocity kicks caused by the RF field near the 
extraction electrode, which lead to different focal lengths 
of the Einzel lens (chromatic abberation). Therefore in 
a future setup the ion lens will be designed to minimize 
chromatic abberation. If in our simulation the RF field 
is switched off before the ion has reached the end caps, 
e.g. 6 periods after start of the extraction, we observe a 
slightly larger spatial spread of the ions after the trap 
(Fig. [51 solid) and before the ion lens (Fig. |B1 solid) but 
the velocity spread in axial direction is strongly reduced 
from Avjv = 5 x lO^'^ to Avjv = 6 x 10^^ resulting in a 
spot size of around 1 nm at the target plane (Fig. [7|3) . In 
order to guarantee the latter velocity spread the power 
supplies have to be stabilized to a level of about 0.5mV.^ 
If the position of the target plane is now changed within 
1 /im the spot size changes drastically (Fig. Ctl-f)- Be- 
side this focussing effect, the lens has an additional ad- 
vantage: The ions move always towards the lens axis. A 
small displacement and misalignment of the ion trap will 
not affect the impact position at the target and makes 
the system insensitive against expected mechanical vi- 
brations and thermal drifts. 

The simulated variation of lateral positions of the ion 
at the target's location as a function of the ion temper- 
ature is shown in Fig. |S1 A displacement below 1 nm 
can be expected for an ion with a temperature below 
0.2 mK if spherical abberation at the ion lens can be 
avoided. This is routinely achieved with laser cooling 
techniques. The simulations also show that a preadjust- 
ment between the trap and the lens within a few /im 
is necessary to achieve a spot size in the order of a few 



* This is quite challenging especially as the voltage has to 
be switched. 
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Fig. 5 Trajectories of 10 ions during extraction with RF 
field switched on (dashed curve) and switched off 6 periods 
after start of extraction (solid curve), respectively. The ion 
trajectory starts in the center of the trap and is directed to 
the upper part of the figure. The end cap electrode of the 
Paul trap is located at 5 mm axial position. Note that the 
radial axis is scaled in nm whereas the axis in the extraction 
direction is scaled in mm. 
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Fig. 6 Trajectories of 10 ions at a distance of 90 mm from 
the trap center (dashed and solid curves as in Fig. |^ . Axial 
position zero denotes the place of the center electrode of the 
Einzel lens. The ions follow linear trajectories from the trap 
and are directed to the upper part of the figure where the 
focus is located at ~ 300 fim axial position. 



Fig. 7 Spot diagram of 500 ray traced ions at the target 
plane. Initial temperature of the single ion inside the trap is 
0.1 mK. The ion is focused with an Einzel lens (located at 
z~ 90 mm) down to a spot size of 20 nm (rms) if the RF 
field is constantly applied and 0.4 nm (rms) if switched off 
6 periods after start of extraction, respectively: a) RF on, 
0.5 fim before focal plane, b) RF on, at focal plane, c) RF 
on, 0.5 after focal plane, d) RF off, 0.5 fim before focal 
plane, e) RF off, at focal plane, f) RF off, 0.5 fim after focal 
plane. 

nm. The focal depth of the lens is about 1 fim and re- 
quires a preadjustnient of the lens voltage. This can be 
done offline using a low energy ion gun combined with a 
detection system, see Section |5| 

For ion beam optics at very low kinetic ion energies 
AC stray magnetic fields as produced by electric engines 
or radio waves are very serious. Fortunately, the ion in- 
side the trap is shielded due to the RF trap field. Sim- 
ulations show that AC stray magnetic fields have to be 
reduced below 0.5 f/T in order to reach nm resolution. 

5 AFM set-up with integrated hollow tip 

With a focal depth of the Einzel lens of about 1 /im 
the distance to the target needs accurate preadjustment. 
Additionally, lateral positioning of the doping sites with 
respect to surface structures, e.g. markers or electric 
gate assemblies is required. The alignment accuracy of 
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piezo resistive cantilevers allows imaging of the 

surfaces' topography in conventional contact mode with 
a resolution in z better than 1 nm and also in non-contact 
mode employing integrated thermal bimorph actuators 
|32|. The Einzel lens is combined with the AFM can- 
tilever mount (Fig. to achieve a high correlation and 
stability between AFM tip and focus position of the lens. 
In contrast to T the AFM tip has no direct ion optical 
function any more and will not be hit by the beam. To 
avoid any kind of contamination during the loading of 
the ion trap, an aperture size in the range of 30-50 nm 
is useful. 



0.1-1 : 1 : 1 : , : , . p- 

200 400 600 800 1000 
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Fig. 8 Calculated rms spot size at the sample as a function 
of ion temperature, RF drive kept on (•) and switched off 6 
periods after start of extraction (o), respectively. Each data 
point is based on 500 simulation runs. 




Fig. 9 Pierced AFM-tip consisting of SiN. The hole in the 
hollow pyramid is fabricated by focused ion beam (FIB) 
milling. The crater shaped pit on top of the pyramid ends 
in a 30 nm hole. 

a few nm is desirable for functional device realization. 
The only known systems that allow addressing with sub- 
nm resolution are atomic-force microscopes (AFM) or 
scanning-tunneling microscopes (STM). 

We take advantage of a recent breakthrough in AFM 
technology: The cantilever probe consists of a SiN-tip 
with a hollow nitride pyramid with a drilled hole (<30 nm 
diameter) fabricated by focused ion beam (FIB) milling 
(Fig. ini) 1221 ■ For this process the Rossendorf IMSA- 
Orsay Physics FIB was applied at an energy of 25 keV 
and a spot size of about 25 nm. A faraday cup under 
the tip was installed as an end point detector to stop 
the milling just in time as well as to image the hole by 
the transmission current signal. The diameter of the en- 
trance of the milled hole was about 140 nm and the exit 
was in the range of 30 nm. The AFM operation with 



6 Application of single ion implantation for 
quantum computing 

For the realization of a solid state quantum computer 
several technical issues have to be solved. The major re- 
quirement is the existence of well defined qubits which 
can be reliably fabricated and controlled. Especially for 
solid-state systems, proper choice of the qubit is a chal- 
lenging task because of decoherence caused by strong in- 
teractions within the bulk material. The most important 
issue at the moment is to proof the scalability of those 
qubit systems. There are several possible approaches like 
qubits based on quantum dots, superconducting nanos- 
tructures and nuclear spins in semiconductors [T1l33ll34j . 
The latter were envisaged as promising candidates be- 
cause of very long coherence times (up to seconds) . 
However, the readout with conventional NMR technol- 
ogy gives only access to ensemble spin states and thus 
causes severe scalability problems IHEl- Recently, it was 
shown, that the single electron spin of a NV center in di- 
amond can be read out optically. As sketched in Figure 
lb, the electronic level scheme exhibits a strong dipole- 
allowed E ^ A optical transition with a zero-phonon 
line at 1.945 eV. A second essential feature of the NV- 
defect is a paramagnetic ground electronic state which 
originates from two unpaired electron spins |37II38II39| . 
The dipolar interaction between these electronic spins 
leads to a zero-field splitting between rus = and rris = 1 
ground state spin sublevels. The readout of the electron 
spin is based on the fact that under excitation of the op- 
tical dipole transition (Fig. ^) the amount of scattered 
photons is sensitive to the internal spin state |40| . 

A two qubit quantum gate was demonstrated us- 
ing the hyperfine coupling between nuclear and elec- 
tron spins associated with a single defect Scaling 
up to a multitude of defects is impeded because of the 
local nature of the hyperfine interaction 41^ . Hence long 
range dipolar coupling between spins is more appropri- 
ate for scalable QC. Note that only neighboring qubits 
will be coupled substantially due to the dependence 
of the coupling. Single qubits may be encoded in elec- 
tron spin states of single NV defects and fast quantum 
gate operations will be performed using electron spin 
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dipole-dipole coupling between closely spaced defects.^ 
As opposed to the hyperfine coupling, now, distances 
of 1 nm are sufficient for achieving a resolvable cou- 
pling of ]7Rabi/27r — 50 MHz. Such interqubit spacings 
are technically feasible by the above described method. 
CNOT and other quantum gate operations between ad- 
jacent electron spins may be implemented using ESR 
analogues of well developed NMR pulse sequences 42, 
I43| . As an example, refocusing techniques using broad 
pulses resonant to all spins may be used for switching 
off all interactions between any spins. In order to selec- 
tively enable the interaction between specific adjacent 
spins an additional spectrally narrow pulse sequence has 
to be inserted between the refocusing pulses. The coher- 
ence time of spins associated with defects in diamond is 
amazingly long with a few tens of microseconds at room 
temperature The two qubit gate operation timescale 
is limited by the coupling strength in the range of tens 
of MHz. This will allow to reach the quantum error cor- 
rection threshold of TgatejT^ — 10"^. 

7 Conclusions 

We have presented a novel method for deterministically 
implanting single ions with sub-nm resolution. Key ele- 
ments of our method are a linear Paul trap serving as 
an ion source in combination with sophisticated laser 
cooling techniques such as EIT cooling to the vibra- 
tional ground state. This atom optical element consti- 
tutes an ideal deterministic point source of ions with 
the ultimate brightness only limited by the Heisenberg 
uncertainty principle. Our simulations predict a lateral 
positioning accuracy of the ions of less than 1 nm for 
an ion-temperature of 0.2 mK reached by laser cooling 
techniques. 

In the future we will optimize the design of the ion 
trap. The ion optical system has to be designed in such 
a way that chromatic abberation is minimized. Exper- 
iments will be performed with various doping ions in 
order to explore a broad spectrum of applications. The 
application of our technique can also be extended to 
molecular or cluster particles. 
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